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ABSTRACT 



We have constructed a new set of isochrones, called the Y'^ Isochrones, that 
represent an update of the Revised Yale Isochrones (RYI), using improved opacities 
and equations of state. Helium diffusion and convective core overshoot have also been 
taken into consideration. This first set of isochrones is for the scaled solar mixture. 
A subsequent paper will consider the effects of a-element enhancement, believed to 
be relevant in many stellar systems. Two additionally significant features of these 
isochrones are that (1) the stellar models start their evolution from the pre- main 
sequence birthline instead of from the zero-age main sequence, and (2) the color 
transformation has been performed using both the latest table of Lejeune et al., and 
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the older, but now modified, Green et al. table. The isochrones have performed well 

under the tests conducted thus far. The reduction in the age of the Galactic globular 
clusters caused by this update in stellar models alone is approximately 15% relative to 
RYI-based studies. When the suggested modification for the a-element enhancement 
is made as well, the total age reduction becomes approximately 20%. When post- RGB 
evolutionary stages are included, we find that the ages of globular clusters derived 
from integrated colors are consistent with the isochrone fitting ages. 

Subject headings: globular:clusters:general - stars:abundances - stars:evolution - 
stars: interiors - cosmology 

1. Introduction 

Isochrones are defined as the locus of coeval (equal age) points on the evolutionary tracks of 
stars of different masses in the Hertzsprung-Russell Diagram (HRD). An older isochrone has a 
fainter and redder main-sequence turn-off (MSTO), because the brighter and bluer massive stars 
evolve and die earlier. For the past four decades, astronomers have used this basic fact in various 
ways to derive the ages of star clusters and galaxies. 

The first systematic application of the isochrone method was made to NGC 188 MSTO by 
Demarque and Larson (1964), as noted by Sandage and Eggen (1969). The widespread adoption 
of the isochrone technique has been tremendously helpful in understanding the formation and 
evolution of the Milky Way and its components, and the chronological sequence of astrophysical 
processes therein (see VandenBerg, Bolte, Sz Stetson 1996; Sarajedini, Chaboyer, & Demarque 
1997). Today, isochrones are the most powerful means of measuring the ages of star clusters, and 
thanks to numerous recent improvements in the input physics, we may now claim to know the age 
of the Milky Way better than ever before. This purely theoretical achievement helps setting one 
of the most important constraints on cosmology, because the age of the universe is a distinctive 
product of any cosmological model. 

A second-generation isochrone technique extends beyond the MSTO up to the tip of the red 
giant branch (RGB), to match more stars in the HRD (Iben 1974; Demarque & McClure 1977). 
The use of the extended isochrone, now standard, confers at least two significant advantages over 
the use of only the MSTO. Firstly, isochrone fitting is no longer in principle sensitive to the 
uncertainty in distance measurement. This advantage, however, has not practically been much 
appreciated because the uncertainties in the stellar atmosphere models and in the convection 
approximations still make it difficult to reproduce accurately the colors of red giants. Secondly, 
it allows us to build realistic population models for clusters and galaxies from their initial mass 
functions. This technique, called "evolutionary population synthesis" (EPS), was pioneered by 
Tinsley (1980). 
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In spite of the attendant difficulties, attempts to reproduce the colors and spectra of clusters 
and galaxies through the EPS technique are increasingly successful, and especially so when 
empirical stellar spectral libraries are resorted to (eg. Larson and Tinsley 1978; Gunn, Stryker, 
&: Tinsley 1981; Bruzual 1983; Pickles 1985). The usage of isochrones has increased steadily. 
Notably, the recent studies of Spinrad et al. (1997) and Yi et al. (2000) have demonstrated that 
precision-isochrones can be used to constrain the galaxy formation epoch directly by estimating 
the ages of high-redshift galaxies in their youth. Such studies have demonstrated that population 
synthesis studies with precision-isochrones may even differentiate between various cosmological 
models. 

Isochrones need to be updated when significant improvements are made in input physics. 
The Yale group (eg. Green, Demarque, & King 1987) is one among several involved in such 
activity. Other notable recent updates include Schaller et al. (1992), Straniero, Chieffi, & Limongi 
(1997), Girardi et al. (2000), and VandenBerg et al. (2000). The Yale group's prior pubhshed 
set of isochrones, the Revised Yale Isochrones ( preen et al. 1987 : hereafter RYI), are out-of-date. 



A number of specialized studies have been published over the years by the Yale group which 
have updated the physics of stellar models (Guenther et al. 1992, 1996; Chaboyer & Kim 1995; 
Guenther & Demarque 1997; Yi et al. 1997; Chaboyer et al. 1998; Heasley et al. 2000). We 
note especially the importance of substantial improvements in the opacities and equation of state 
(mostly originating from the OPAL group: [glesias fc Rogers 19961) . Because of the need for a 



comprehensive set of isochrones for research in stellar populations, we hereby present a grid of 
Yale isochrones using up-to-date input physics and parameters. 

Our new isochrone set ("the Isochrones", after the Yonsei-Yale collaboration) covers a 
wide range in metallicity and age. A wide Z-coverage in particular is useful in constructing EPS 
models for complex populations, such as elliptical galaxies. The ages in the Isochrones are 
computed starting from the pre-MS birthline rather than from the "zero-age" MS (ZAMS). This 
allows us to build realistic young isochrones that are mainly defined by pre-MS stellar models. In 
this respect, they should be useful in population studies of youthful systems, such as young open 
clusters. 

In this paper (the Y^ isochrones-solar mixture), all isochrones are based on a chemical 
composition that has been scaled to the solar mixture. The case of a-element enhancement will 
be considered in a subsequent paper. This is because we believe that we do not as yet have a clear 
understanding of it as a function of mass and of metallicity, particularly in the high metallicity 
regime. If a-element enhanced isochrones is desired, we recommend the use of a correction formula 
similar to the one suggested by Salaris et al. (1993) with this set of isochrones. This correction 
formula has proven reliable in low metallicity systems, such as the globular clusters (Chaboyer et 
al. 1992), when compared to models actually constructed with opacities for a-element enhanced 
mixtures. Subsequent work will present Y^ isochrones that include a-element enhancement based 
on more up-to-date opacities. Users will then be able to construct isochrones for the desired values 
of the a-element enhancement via simple interpolation. The corresponding post-RGB models will 
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also be available presently. 



2. Input Physics and Parameters 

The physics used in this study for the construction of stellar models is summarized in Table 1 
and discussed below. 



2.1. Microscopic Physics 

The OPAL opacities are the most widely used Rosseland mean opacities today. Lately, the 
group has released newer tables for the Grevesse and Noels (1993) solar mixture that include the 
effects of seven additional heavy elements. These new calculations also include the effect of some 
changes in physics ( Rogers fc: Iglesias 1995 ; [glesias fc Rogers 1996| ). All these effects have been 



taken into account. We use the OPAL opacities for the temperature range of logT > 4.1. For 
logT < 3.75, we use the Alexander & Ferguson opacities (1994). In the temperature region of 
3.75 < logT < 4.1, we use a value linearly interpolated between these two sets of tables. The 
conductive opacity is always included when logT > 4.2 and log/? > 2 x logT — 13, where p is 
density. We adopt the work of Hubbard & Lampe (1969) for logp < 6.0 and the work of Canuto 
(1970) in the relativistic regime where log p > 6.0. Deep inside the star where the OPAL opacity 
table fails to cover, the electron scattering is the dominant mechanism. 

The equation of state was taken from Rogers, Swenson, & Iglesias (1996), i.e., OPAL EOS. 
Beyond the boundaries of the table, the standard Yale implementation with the Debye-Hiickel 
correction was used ( puenther et al. 1992| ; [Chaboyer &: Kim 1995 ). The most popular equations 



of state have been simple Saha solvers, which we used when logT < 3.7. As the importance of the 
Coulomb interaction has been acknowledged, the effect has been included using the Debye-Hiickel 
approximation. Care has to be taken, however, that the Debye-Hiickel correction is only applied 
when A < 0.2 (see page 23-26 for definition of A and discussion in Rogers (1994)). For masses 
of O.75M0 and above, A is always less than 0.2 in the post-ZAMS stellar evolution models. 
Thus, models constructed using the OPAL equations of state and those of the usual Saha solver 



together with the Debye-Hiickel approximation are similar ( Chaboyer fc Kim 1995 ). Note that 
the treatment of Coulomb interaction in the OPAL equations of state are valid for values of A up 
to 3, more than an order of magnitude better than the simple Debye-Hiickel formula in common 
use. At very high temperature and/or pressure, where the OPAL EOS tables do not cover, we 
fall back to the YREC built-in EOS, which follows Cox & Giuli (1968), and is described in detail 
in Appendix A of Prather (1976). Here the pressure ionization is the dominant effect, i.e. fully 
ionized medium. Between tables, we have interpolated linearly. 

To include helium diffusion in the calculation, we have employed Loeb's formula (Bahcall, 
& Loeb 1990; Thoul, Bahcall, & Loeb 1994). For a more detailed discussion of diffusion, see 
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Chaboyer et al. (1992)| where the work of [Thoul et al. (1994)| and [Michaud and Proffitt (1993) 



have been compared. There is an uncertainty in the diffusive coefficient of each metal element. 
Thus, common practice assumes that all heavy elements diffuse with the same velocity as fully 
ionized iron. Guenther &: Demarque (1997) present the results of such computations. This study 
includes only helium diffusion. 

The energy generation routines of Bahcall & Pinsonneault (1992) and the cross sections 
listed in Bahcall (1989) have been used. In addition, the cross sections for the pp reactions, the 
-Be''-proton capture reaction, and the hep reaction have been updated (Bahcall and Pinsonneault, 
priv. comm.). These routines, kindly provided by Bahcall (priv. comm.), also include coefficients 
for weak, intermediate and strong electron screening, based on the Salpeter (1954) theory of weak 
screening (which is valid in most cases in stellar evolution calculations), and its extension to 
intermediate and strong screening (DeWitt et al. 1973; Graboske et al. 1973). Neutrino losses 
from photo, pair and plasma neutrinos are included in the energy generation following the work of 
Itoh et al. (1989). 



2.2. Solar Calibration 

Our stellar models have been calibrated against the sun. These models are based on the 
solar mixture of Grevesse & Noels (1993). Grevesse &: Noels made a major effort to determine 
more accurate CNO & Fe abundances, and now the solar Fe abundance agrees with the meteoritic 
abundance better than ever before. Their new solar metal-to-hydrogen ratio is {Z/X)q = 0.0245, 
while the previous value from Anders & Grevesse (1989), based on the meteoritic Fe, was 0.0267. 
The value we have finally adopted is 0.244 from the even more up-to-date value from Grevesse, 
Noels, & Sauval (1996). 

The model that best matches the solar properties {Lq = 3.8515E33 ergss"^ and Rq = 
6.9551E10 cm) at the generally accepted age of the sun (4.50 Gyr) is of {X,Z)o = (0.7149,0.0181), 
i.e., {Z/X)q = 0.0243. Our solar model achieves the following agreements: T^s/Tcsq = 0.999936, 
log (L/Lq) = 0.0000000185, and log (R/Rq) = -0.0000006701. At the solar age, the surface 
composition was (X, Z) = (0.7463, 0.0181). The mixing length parameter of l/Hp = 1.7432 has 
been used to produce this match and thus for all other models. 

Note that we have used the slightly larger value of Lq (3.8515E33 ergs s~^) rather than the 
one listed in Livingston (2000), i.e., 3.845E33 ergs s~^. This 0.17% difference in luminosity is due 
to the uncertainty in the measured value of the solar irradiance. Guenther & Demarque (1992) 
state that "the luminosity has been determined from solar constant measurements from space 
on both the Nimbus 7 and the SMM satellites (Hickey & Alton 1983). ERB-Nimbus measures 
1371.6 ± 0.765Wm2 and SMM/ACRIM measures 1367.7 ± 0.802Wm2 which yields luminosities 
of 3.846E33 ergs~^ and 3.857E33 ergs~^, respectively. As we cannot expertly base a preference, 
we merely take the average of the two and adopt the solar luminosity in all our models to be 
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Lq = 3.8515E33 erg s~ ." The accuracy of the Nimbus 7 measurements has more recently been 
estimated to be 0.5% (Lee et al. 1995). We simply adopt the value of Lq from Guenther &; 
Demarque (1992) in this study. 

2.3. Chemical Abundance and AY/AZ 

We have set the initial chemical composition to (Y,Z)q = (0.23,0.00). Our solar calibration 
(described above) suggests the initial solar chemical composition of {Y, Z)q = (0.267025, 0.018100). 
This indicates AY/AZ ^ 2.0, i.e., 

y = 0.23 + 2Z. (1) 

This slope is the same as that from VandenBerg et al. (2000) and only slightly smaller than that 
(2.25) of Girardi et al. (2000). Given the inevitable differences in codes and input parameters, 
we regard this as a good agreement. Table 2 displays the chemical compositions selected for this 
study. 

It should be noted that AY/AZ is not a precisely determined quantity. In this study, it has 
been determined based on two values: the initial composition and the solar composition. Various 
techniques have indicated differing values, but most estimates cluster between 2 and 5 (sec Pagel 
Sz Portinari 1998 for review). Adoption of a slightly different value does not greatly modify the 
stellar evolution when Z is low. However, it can lead to unrealistic stellar models when Z is 
very large (> 2Zq), because helium abundance is the prime factor in setting the pace of stellar 
evolution. In this sense, our extremely metal-rich models based on a crude value of AY/AZ may 
not have the same accuracy as their metal-poor counterparts. 

2.4. Correction for cc-Element Enhancement 

The isochrones presented in this paper are intended for solar-type populations that do not 
show signs of a enhancement. For populations with strong signs of a enhancement, such as 
metal-poor globular clusters, we recommend the correction formula similar to the one provided by 
Salaris, Chieffi, & Straniero (1993): 

Z = Zo(0.638/a + 0.362), (2) 

or, one more appropriate to the Grevesse Sz Noels (1993) mixture, which is: 

Z = Zo(0.694/c,-F 0.306), (3) 

where fa is the chosen a enhancement factor, Zo is the actual metallicity of the target, and Z 
is the metallicity of the appropriate non-a-enhanced isochrone. The modification of the original 
Salaris et al. formula is necessary because of the change in the mixture adopted. This effect is 
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probably small, as indicated by Salaris et al. (1993). But we also emphasize that such corrections 
are likely to be invalid at the higher metallicities. 

Upon availability of our second set of isochrones (with a enhancement included), we will 
validate this issue and users will be able to generate the appropriate isochrones for any reasonable 
a enhancement factor through a simple interpolation. 



2.5. Convective Core Overshoot 

Convective core overshoot (OS), the importance of which was first pointed out by Shaviv & 
Salpeter (1973), is the inertia-induced penetrative motion of convective cells, reaching beyond the 
convective core as defined by the classic Schwarzschild criterion. Stars develop convective cores if 
their masses are larger than approximately 1-2 Mq, typical for the MSTO stars in 1 - 5 Gyr-old 
populations, depending on their chemical composition. Since the advent of the OPAL opacities, 
various studies have suggested a modest amount of OS; that is, OS 0.2 Hp for clusters of age 1 
- 2 Gyr but OS ~ 0.0 - 0.1 for older clusters (4-6 Gyr), where is the pressure scale height 
(Stothers 1991; Demarque, Sarajedini, & Guo 1994; pinescu et al. 1995 ; Kozhurina-Platais et al 



1997) 



In this study, we have adopted OS = 0.2 for younger isochrones (< 2 Gyr) and OS = 0.0 
for older ones (> 3 Gyr). This is based on the observational studies listed above. It should be 
noted, however, that such tests of the models were performed mainly against the approximately 
solar abundance populations, and thus our adoption of the OS parameter may not be valid in 
chemically different populations. 

For young isochrones with OS, we first found the critical mass above which stars 

continue to have a substantial convective core even after the pre-MS phase is ended, by inspecting 
the stellar models with no OS. The mass interval between our adjacent models being 0.1 Mq, 
the actual W^^f;" is between [M'Sit listed) and (M^^;;^'^ - 0.1). Then, we constructed heavier 
stellar models with OS included and used them in the construction of young isochrones. For old 
isochrones, we used stellar models with no OS regardless of mass. OS has many effects on stellar 
evolution (see Stothers 1991). Among the most notable are its effects on the shape of the MSTO, 
on the rate of evolution near the MSTO and the subgiant phase, and on the ratio of total lifetimes 
spent in the core hydrogen burning stage (MS) and in the shell hydrogen burning stage (RGB). 
The impacts of such effects on the isochrone and on the integrated spectrum have been discussed 
by Yi et al. (2000). 



3. Construction of Stellar Models 
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3.1. Evolutionary Tracks 

We have used the Yale Stellar Evolution Code (YREC) to construct these stellar models. In 
order to construct the isochrones for the age range 0.1-20 Gyr that extend to the tip of the RGB, 
we have constructed evolutionary tracks for mass 0.4 - 5.2 Mq, from the pre-MS birthline to the 
onset of helium burning. In addition to these full isochrones, we also provide younger isochrones 
that reach only up to the upper MS, as explained in §5. 

3.2. Initial Models at Pre-MS Birthline 

An important improvement of the isochrones over RYI is the fact that these stellar models 
begin, not at the classic ZAMS, but at the deuterium MS, also called the "stellar birthline", where 
stars initially become visible objects (Palla and Stahler 1993), during the pre-MS stage. Thus, 
these models include pre-MS evolution. 

The pre-MS lifetime, denoting the time taken by a star to evolve from the birthline to the 
ZAMS, is a strong function of stellar mass, varying from less than 1 Myr to about 200 Myr for stars 
in our mass range. It is about 43 Myr long for a solar mass model. This means that the ZAMS, 
as hitherto defined, has an inherent age spread of about 200 Myr, making it simply inappropriate 
to use models starting from the ZAMS to match the HRDs of young clusters. Because of this, 
isochrones that label ZAMS stars as "zero-age" are bound to mismatch the lower part of the 
observed HRDs of young clusters that contain pre-MS stars (e.g. Patten & Pavlovsky 1999). This 
effect is negligible in old (> 1 Gyr) populations, causing age underestimates of only 1% for 10 Gyr 
populations, but it can be significant in the analysis of young open clusters. We therefore suggest 
that ages be henceforth defined relative to the birthline rather than the ZAMS. 

In view of the expanding utility of isochrones beyond the treatment of old systems, in 
recognition of the increasing importance of very young systems, and to provide a unified treatment 
for open, globular and perhaps other star systems, we have decided to include pre-MS evolution 
in this set of isochrones. 

This raises the issue of where and how the models should be started off. Is there some 
equivalent of the ZAMS on the Hayashi track? One cannot begin the evolution at an arbitrary 
point on the pre-MS, since the radius and the central density of stars, both important quantities, 
especially from the viewpoint of stellar rotation, but also in a more general sense, change rapidly 
during this phase. This deprives us of the pause on the MS, when hydrogen burning begins, that 
serves as a natural starting point for models relevant to old stellar populations. An equivalent of 
this pause is provided by the deuterium MS, or stellar birthline, as defined by Palla and Stahler 
(1991)0. Its use for these purposes seems both apposite and physically meaningful. For the 



^We note that, in contrast to the situation with the MS, this is a theoretical construct, and there is not, at least 
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low-mass stars in this study, the birthUne is distinct from and above the MS, located just above 
the T Tauri stars in the HRD. 

We have constructed a series of birthline models for low-mass stars on the birthline. These 
models are then evolved forward in time along the pre-ms tracks, onto the MS and beyond, ending 
at the RGB tip, but thus making available isochrones intended for quite young systems as well, 
perhaps valid even for 10 Myr-old systems. 

Operationally, these models were started off from polytropic models that were placed even 
higher than the birthline in the HRD. However, polytropic models are quite different from real 
stellar models in their thermal structure. Thus they need to be relaxed until their thermal 
structures are correct. This is accomplished in YREC and the resulting structures are evolved 
down their Hayashi tracks until they satisfy the mass-radius relationship of Palla & Stahler (1991), 
at which point they are assigned age zero (see also Barnes & Sofia [1996]). This pre-birthline 
evolutionary phase is a convenient way of calculating physically self consistent initial models, and 
has no special physical significance since it does not include deuterium burning in the definition 
of the birthline. Deuterium burning is however implicitly included in the post-birthline evolution 
in the reaction network of YREC. The initial deuterium gets destroyed, but the contribution of 
deuterium burning to the total energy production is negligible compared to the contribution from 
gravitational contraction. Thus, the starting models have both the correct, or at least appropriate, 
radius and also internal structure^. The chemical composition is now checked, and adjusted to 
be pristine. Models with differing chemical compositions may be constructed by rescaling the 
corresponding solar mixture models. A series of these models have been constructed for different 
masses. These models are then evolved onto the MS and beyond. 



4. Color Transformations 

Theoretical properties ([Fe/H], L, T^s) have been transformed into colors and magnitudes 
using the color transformation tables of Lejeune, Cuisinier, & Buser (1998, hereafter LCB) and of 
Green et al. (1987, hereafter GDK), which was used in the RYI. Both tables are semi-empirical: 
the LCB table is based on the latest Kurucz spectral library (1992) and the GDK table is based 
on the older Kurucz library (1979). Both tables have been substantially modified to match the 
empirical stellar data better than their forebears (Kurucz tables) and extended to low temperatures 
based on empirical data and low temperature stellar models. For the present calculations, we have 
used the updated BaSeL-2.2 version of the LCB library^, for which the calibration of the cool giant 



not yet, an observational definition of tfie deuterium MS. 

''Tfiese issues are particularly important in studies of the early rotational evolution of stars because of the 
importance of the stellar moment of inertia. 



The BaSeL-2.2 tables are available only electronically at ftp://tangerine.astro.mat.uc.pt/pub/BaSeL/ 
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models in the parameter ranges 4500 K < T^s < 6000 K and 2.5 < log^ < 3.5 has been revised. 

In the process of comparing isochrones to observation, the adoption of an up-to-date color 
transformation table is as important as updating input physics in the construction of stellar 
models. The difference between the LCB table and the GDK table is substantial, as shown in 
Figures 1 and 2. Figure 1 shows the comparison of bolometric correction (BC) for logg = 4.5 MS 
stars. The overall agreement is good. At solar metallicity and temperature (vertical line in the 
top right panel), BC from the LCB table is somewhat higher than that from the GDK table. For 
example, BCq is —0.08 and —0.109 in the GDK table and in the LCB table, respectively. We have 
normalized these two tables to the Sun so that the visual absolute magnitude of the Sun {Mvq) 
becomes 4.82 (Livingston 2000), regardless of the table used. Then, the following definition of BC: 

BC = MboiQ - Mv (4) 

yields M^oiq = 4.711 and 4.74 for the LCB table and the GDK table, respectively. Because My is 
defined as: 

My = MboiQ - 2.51ogL/L0 - BC, (5) 

the two tables yield the same V magnitudes. 

Figure 2 shows how seriously the two color tables difi"er from each other. The original GDK 
table reaches only down to the effective temperature of 2800K and up to 20000K. In order to 
cover the tip of the metal-rich RGB {Z > 0.06), we have extrapolated this table down to 2500K 
when necessary. A few of the youngest and most metal-poor isochrones reach above effective 
temperature 20000K. In that case, we have used the colors from the LCB table. This crude remedy 
did not cause any noticeable discontinuity in the isochrones. 

The LCB table is obviously more up-to-date than the GDK table. Yet, it is not always 
our preference. In some circumstances, the isochrones with the GDK table match the empirical 
data better. However, without knowing the true parameters (metallicity, reddening, and distance 
modulus), it is difficult to choose one over the other. Thus, we have decided to provide our 
isochrones in two formats based on both tables. Users are recommended to try both. The color 
transformation is for the filter systems of {UBV)johnson{BI) Cousins in the GDK table and of 
{UBV)johnson{RI)cousins{JHKLM)ESO in the LCB table. 

5. Results 

5.1. Stellar Evolutionary Tracks 

Figure 3 shows a sample set of stellar evolutionary tracks for the solar composition and for 
the mass range of 0.4 - 5.0 Mq. The line connecting the diamond symbols is the birthline for 
this composition. The first position of each track in the HRD is somewhat uncertain because it 
depends on the choice of the first time step in the numerical computation. This uncertainty can 
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be alleviated if we force the first time step to be very small. However, we have decided not to pay 

too much attention to such details in the first steps of the pre-MS stage because we believe that 
other uncertainties in the position of the birthline and in observational data are much larger. In 
this sense, our youngest isochrones below lOMyr are less certain than their older counterparts. 

5.2. Isochrones 

Figure 4 shows a set of isochrones for Z = 0.02 and ages of 1 Myr through 20 Gyr. The young 
isochrones (1, 2, 4, 8, 10, 20, 40, 60, 80 Myr) are complete only to the MS and thus focused on the 
pre-MS stage. As mentioned earlier, one of the most notable features of the Y"^ Isochrones is the 
inclusion of pre-MS stellar evolution. 

The information that we provide in the Y"^ Isochrones is listed in Table 3. The last three 
columns are luminosity functions (LFs); namely, the number of stars in the box defined by the 

mass grid given assuming there were initially 1000 stars in the mass range of 0.5 - 1.0 in Mq, 
following the concept introduced in the RYI. The bottom right panel in Figure 4 shows the 
evolution of LF in the case of IMF slope x = 1.35 (Salpeter index). 

Figure 5 compares the Isochrones and RYI. One can see the substantial change in the 
stellar models in the top left panel. The other three panels show the isochrones in the observer's 
HRD. The effects of the use of the LCB color transformation table are clear. With a close 
inspection of the lower MS of some isochrones based on the LCB table, one will find some 
discontinuities m.U — B and B — V colors (e.g., lower left panel). These are partly real because 
they are also present in the empirical data, although with a smaller amplitude. The discontinuities 
in the isochrones get amplified by the fact that the empirical data are scarce especially in terms of 
metallicity and gravity grids. Besides, the adjustment (correction) between the Kurucz spectral 
library (the foundation of our color tables) and the existing low temperature libraries made by 
LCB is not as feasible as we had hoped for because of the large and apparent difference between 
them (At times, the flux difference is as large as an order of magnitude). 

Similarly, peculiarly curved features on the RGB tip are shown in the LCB-based isochrones 
(inU — B and B — V colors) . These features are caused by the non-monotonic temperature-color 

relations that are present in the empirical data. Lejeune et al. (1997, Figures 6 & 14) illustrate 
these effects. Whether the empirical stellar data that show such non-monotonic temperature-color 
relations are being correctly interpreted and implemented in the color transformation table 
remains to be determined through more rigorous data sampling. Because of the limitation of the 
adopted calibration procedure by LCB, the corrected LCB colors still differ from the empirical 
data by as much as 0.1 mag 'm.U — B and B — V (see for instance Figure 14 in LCB). In this sense, 
our LCB-based isochrones are subject to such uncertainties. Readers are referred to Lejeune et al. 
(1997; 1998) for more information. Despite the magnitude of the disagreement in the colors in 
this temperature range (often up to Imag!), this may not be an important issue in the study of 



- 12 - 



most Milky Way globular clusters, because their metallicities are too low to be affected by this 
uncertainty. But it is a serious problem in metal-rich stellar populations of the kind found in the 
bulges of spiral and elliptical galaxies. 

Part of a sample isochrone is shown in Table 4. The tables in their entire form are available 
electronically from this paper. Minor corrections may occur to the tables after publication. The 
updated versions will be available from the authors upon request or directly from our WEB 
site0. We also provide a FORTRAN interpolation routine that works for metallicity and age 
interpolation via the same WEB site. 



5.3. Bump on the Red Giant Branch 

The evolution in luminosity on the RGB is not always a monotonic function of time. When the 
depth of the convective envelope increases, some processed material is distributed over the entire 
convective region, and CNO processed isotopes are mixed to the surface. Stars are rejuvenated by 
such mixing because unprocessed hydrogen is mixed from the envelope into the interior. The first, 
and largest, of these partial mixing phenomena on the RGB, is called "first dredge up", an event 
which takes place when the star's convection zone reaches its maximum depth soon after it reaches 
the giant branch. The convection zone depth then decreases (in terms of mass fraction) as the 
star evolves up the giant branch, leaving behind a composition discontinuity. At some later point, 
the hydrogen burning shell passes through the composition discontinuity. This event gives rise to 
a temporary decrease in luminosity, as shown in Figure 6. This was predicted by Thomas (1967). 
This results in a peak in the luminosity function on the RGB, or "bump", which is observed in 
some globular clusters (see e.g. King et al. 1985; Fusi Pecci et al. 1990). At luminosities higher 
than the bump on the giant branch, stars will have more fuel to burn in the upper part of the 
RGB, because the first dredge-up had previously mixed in fresh hydrogen from the envelope. 

Figure 7 and Table 5 display the predicted position of the RGB bump as a function of age 
and metallicity. Except when the metallicity is extremely low, the RGB bump generally gets 
fainter as age and metallicity increase. We have also plotted in Figure 7-(b) the V magnitudes of 
the RGB-bumps in some globular clusters. The V brightness data are from Table 4 in Fusi Pecci 
et al., while the distance moduli and the metallicities are from Harris (1996). In the case of 47Tuc 
(the second most metal-rich data point), we have adopted a distance modulus of 13.47 from our 
CMD fit described in Figure 13. The most metal-rich data point {V{bump) ~ 17.1) has been 
derived from the NGC6553 data from Zoccah et al. (2001). Adopting [Fe/H]=-0.17 (see §7.3), 
we derived Mv{bump) ~ 1.3, which is also in good agreement with the prediction. 

The prediction of the position of the RGB bump cannot be precise because of the uncertainties 
in the detailed physics at the base of the convection zone and the dynamics of the first dredge-up. 
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When the RGB bump was first identified in 47Tuc (King et al. 1985) the observed RGB bump 
luminosity then differed from the predicted luminosity in the theoretical tracks, and possible 
causes for this discrepancy were explored, principally convective overshoot below the convection 
zone, and opacity uncertainties. Subsequent improvements in the quality of the observational data 
and in the opacities and other physics input have resulted in a decrease in the discrepancy (Pusi 
Pecci et al. 1990). More recently, some researchers have argued that the discrepancy disappears 
when one takes into account the opacity increase due to alpha-enhancement and that the effect of 
convective overshoot is negligible (Cassisi and Salaris 1997; Ferraro et al. 1999). 

We note at this point that other physical processes can play an equally significant role. The 
presence of mixing other than convective overshoot complicates the issue. There is evidence that 
some mixing below the convection zone takes place in red giants, probably the result of meridional 
circulation due to rotation (Sweigart & Mengel 1979). Near the MS, the envelope overshoot 
does not seem sufficient to explain the observed lithium depletion below the convective envelope 
(Pinsonneault 1994). We also note that helioseismological inversions indicate that there is little 
overshoot below the convection zone in the present sun (Basu, Antia, k, Narasimha 1994). For all 
these reasons, our models were calculated without any overshoot below the convection envelope. 
The effects of including alpha-element enhancement on the RGB bump position will be discussed 
in our next paper. 

5.4. Maximum Brightness of the RGB 

In general, the RYI did not provide accurate positions of RGB tips, because Tintil recently 
isochrones were dominantly used for HRD fitting, a matter in which the RGB tip plays a small, 
if any, role. However, isochrones are now used in a variety of studies. One example is a field of 
evolutionary population synthesis, whose aims inchidc the computation of integrated properties 
(eg. colors and magnitudes) of stellar populations. The position of the RGB tip is important 
in computing near infrared magnitudes, as stars near the RGB tip are good infrared sources. 
Therefore, we have paid extra attention in order to define the RGB tip as consistently as possible. 
It should be noted, however, that evolution near the RGB tip remains somewhat uncertain, and 
will remain so as long as mass loss is poorly understood. 

Another important application of accurately modeled RGBs is to use them as distance 
indicators for metal-rich populations. Figure 8 shows that the brightness in V magnitude on the 
RGB reaches a maximum not at the RGB tip but in the middle of the RGB. Because the number 
density of stars in the middle of the RGB is much higher than that of the RGB tip, this gives us 
a large statistical advantage. Besides, this part is more certain than near the RGB tip in terms 
of our understanding of stellar evolution. This effect, however, occurs only when the metallicity is 
high enough, i.e., Z ^ 0.004, because it is mainly an opacity effect. There are currently a great 
number of observational studies that are based on this concept. New predictions are provided in 
Figure 8 and Table 6. 
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6. Luminosity Functions and Photometric Evolution 

Isochrones are convenient building blocks for the the study of evolutionary population 
synthesis (EPS). EPS studies either use the luminosity functions (LFs) provided in the isochrone 
tables (if provided as is the case with our isochrones) or compute the LFs using the mass 
information that is always provided in the isochrones. Generating isochrones is not a trivial 
matter partly because keeping the desired accuracy in mass through several track interpolations 
is not easy. As Yi et al. (2000) mentioned, the mass difference between the bottom and the tip 
of the RGB in the 12Gyr isochrones is merely an order of 0.7% or approximately 0.006 Mq. If 
the mass interpolation is inaccurate in the 5th decimal place in Mq, the LF will lose its desired 
accuracy. Such errors can be checked by inspecting the integrated photometric properties of a 
sample population, as shown in Figure 9. Our isochrones all show clear and smooth luminosity 
(and color) evolution, except where there are expected departures from the monotonicity in the 
color transformation table. Considering this complexity, it is better to use the isochrones rather 
than stellar tracks in most of the EPS studies. 

7. Tests of Isochrones against Empirical Data 

We are encouraged by the fact that the improved isochrones actually provide a better match 
to the empirical data. In this section, we display some of the tests we have made. We have 
attempted to provide a variety of tests to check the validity of different aspects of the isochrones. 

For example, we try to test our isochrones on young open clusters as well as old globular clusters. 
We demonstrate the level of agreement between the isochrones and the MS stars in the field as well 
as in an open cluster. We also present a test of the color evolution of Galactic globular clusters. 
We are pleased to find that our isochrones combined with their LFs reproduce their integrated 
colors quite well at their generally accepted ages. 

7.1. Population I MS stars 

One of the most basic tests is to examine whether our stellar models match the measured 
properties of nearby MS stars. Figure 10 shows a set of stellar evolutionary tracks for Z = 0.02, 
i.e., nearly solar. Overplotted are the observed MS stars from Gray (1992). The match is good: 
not only do the models match the locus of the observed MS, but the match in mass is also 
good except in the very low mass range (M ^ 0.6 Mq). The masses of the faintest three data 
points (M dwarfs) are 0.6, 0.52, and 0.48 Mq according to Gray (1992), while models suggest 
approximately 0.7, 0.6, 0.58 Mq. This disagreement in mass in the low mass range may be due 
to the uncertainties in the mass determination of cool stars or in the low-temperature stellar 
atmosphere models. 
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7.2. Subdwarf MS Stars 

Figure 11 shows the match between our isochrones and the subdwarf MS stars whose distance 
has been determined through HIPPARCOS observations. The data are from Reid (1998), and 
we have fohowed the same criteria for the data selection as VandenBerg et al. (2000) did for 
their Figure 14: i.e., [Fe/H] < -0.55, My > 4.5, and a^/ir < 0.07. Instead of 14Gyr 
isochrones, we are using 12 Gyr isochrones. However, there is practically no difference in this part 
of the isochrones. Since our isochrones do not include the effects of a-element enhancement, we 
have recalibrated to compute the isochrones in this plot using Equation [2]. We have adopted 
[a/Fe]=+0.4 and +0.2 for [Fe/H]< -1.0 and [Fe/H]=-0.5, respectively, following Wheeler, 
Sneden, Truran (1989) and Carney (1996). The agreement is good and is marginally better with 
the isochrones based on the GDK table. 



7.3. Globular clusters 

Globular clusters are the main targets of isochrone applications. They are the playgrounds 
of stars of the same age and metallicity, so we can test and refine our stellar models. Moreover, 
an accurate estimation of their ages has been appreciated as an important route to a better 
understanding of cosmology. 

Figures 12 and 13 shows fits to the metal-rich cluster 47Tuc and the metal-poor cluster M68, 
using our isochrones based on the GDK table. We have used theoretical values of My (RR) from 
Demarque et al. (2000) to derive their distances. These models are based on the HB evolutionary 
tracks with input physics nearly identical to the current isochrones. They also take into account 
the effects of HB morphology in the calculation of Mv'(RR) for a given [Fe/H]. The theoretical 
values of My (RR) from these models are in a reasonable agreement with RR Lyrae luminosities 
derived from the main sequence fitting of globular clusters based on HIPPARCOS parallaxes for 



field subdwarfs (Reid 1997; Gratton et al. 1997). For purposes of illustration, we have shown 



synthetic HB models in the lower panels (see Lee, Demarque, & Zinn 1994 for description). 

In case of 47Tuc, its red HB is approximately 0.15 mag brighter than RR Lyrae level. The 
matches are good if their ages are approximately 12 and 13 Gyr, respectively. Through the tests 
on a few other globular clusters, we have found a reduction in the age estimates of the order of 
20%, and this is when a-enhancement has been considered using the Safaris et al. formula. 

The primary source of this age reduction is the update of the stellar models. Of several ways 
to measure the ages of globular clusters, the most popular one is probably to match the MSTO 
luminosity using isochrones. The Y'^ isochrones are fainter at the MSTO (bluest point on the MS) 
than the RYI by approximately 16% when Z = 0.0004, a metallicity typical for Galactic globular 
clusters. Figure 14 clearly shows this. A given observed MSTO luminosity indicates substantially 
younger ages when the isochrones are used. Also compared in this figure is the result derived 
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from the isochrones of Girardi et al. (2000), whose input physics is quite compatible to ours. The 
agreement with our models is good. 

Table 7 shows the change of the derived ages caused by the use of the new stellar models alone. 
The estimates are the mean values of the age estimates for various values of MSTO luminosities 
within the range of 5 16 Gyrs. Much of this age reduction is caused by the inclusion of the the 
helium-diffusion and the use of the updated equations of states. Our preliminary investigation 
suggests that, in addition to this, the inclusion of alpha-enhancement causes approximately 
another 5% of age reduction. It should also be noted that the change (reduction) in derived ages 
caused by the improvement in the stellar models alone (no alpha-enhancement effects included) is 
reversed in the metal-rich regime (see Table 7). These issues will be discussed more thoroughly in 
our next paper. 

The next test is on NGC6553, a metal-rich Galactic globular cluster (Figure 15). The faint 
data (dots) are from Zoccali ct al. (2001). They are confirmed members of the cluster. Also 
shown are the bright stars {V < 17: crosses) from Sagar et al. (1999). They have been selected 
to satisfy Erroriy — I) < 0.2. These two data sets have offsets both in V magnitude and in V-I 
by 0.3 mag and 0.15 mag (the Zoccali et al. data being brighter and bluer), respectively. This 
is because the Zoccali et al. data are from the region that suffers less reddening (Zoccali, priv. 
comm.). When we shifted the Sagar et al. data, both sets match one another from the MS to the 
HB clump very well. The magnitudes of the offsets were estimated from the mean positions of the 
HB clump through eye-fit. 

The metallicity of this cluster has been reported to be in the range of [Fe/H]=— 0.5 — hO.l 
(see Sagar et al. 1999), i.e. approximately Z ~ 0.006 - 0.02. We have adopted Z=0.0125 
([Fe/H]=— 0.17) in this study because at this metallicity the LCB-based Y"^ isochrones match 
the upper RGB of this cluster best. The age of the model (9Gyr for this fit, E{y-I)=Q.'d, 
m-M[apparent]=15.78) has been derived from the luminosity difference between the MSTO and 
the mean luminosity of the HB clump. The open diamond shown on top of the HB clump is the 
synthetic HB clump for this age and metallicity. However, the age estimate is sensitive to the 
adopted metallicity; 10 - 8Gyr when Z=0.01 - 0.02. Zoccali et al. (2001) achieve a somewhat 
larger age estimate, 12 Gyr, but this is mainly because they assume a metallicity lower than ours. 
Using LCB-based isochrones, we favor our metallicity and thus a smaller age. 

The match by the LCB-based Isochrones is good. The turn-around in the V magnitude 
would not be as dramatic as observed if the GDK table is used, although GDK-based isochrones 
seem to match the MS better. Considering the importance of this bulge cluster to the study of the 
formation of Milky Way, a more careful analysis, in particular considering a-enhancement (highly 
suspected to be present [B. Barbuy, priv. comm.], is required to give a definite answer. 
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7.4. Open clusters 

If tests against globular clusters teach us about the stellar evolution in the MS and phases 
beyond, those against young open clusters teach us about pre-MS stellar evolution. Several 
studies have noted that the theoretical isochrones are fainter in the lower MS than observed (e.g., 
Subramaniam & Sagar 1999; Patten &; Pavlovsky 1999). This is likely to be evidence of pre-MS 
evolution. 

Figures 16 and 17 show isochrone fits to the young open clusters, Pleiades and IC2391. Both 
LCB- and GDK-bascd isochrones work well for these fittings. The Pleiades data are very well 
represented by the 40-100 Myr isochrones, assuming E{B — V) = 0.04 and m — M(apparent)=5.6 
(Pinsonneault et al. 1998). A good fit requires an apparent distance modulus somewhat larger 
than that derived from the HIPPARCOS data, as was noted earlier by Pinsonneault et al. (1998). 

The overall fit appears to be good, but the data at 5 — F « 0.5 and ^ 1.5 seem to be notably 
brighter (or redder) than the models. If the observational errors are not to be considered large, 
this may indicate the shortcomings of the models. 

The MS data of the younger cluster IC 2391 is matched by the 20 Myr solar-composition 
isochrone. This is reasonably close, if somewhat smaller than the age estimate of Patten & 
Pavlovsky (1999), i.e. 30 Myr. The membership study of this cluster has been performed very 
carefully, and thus the data are of good quality and the sequence is better defined than the data 
shown in Figure 16. The model isochrones are for 4, 10, 20, 40, and 80 Myr of age, while the 
20 Myr isochrone is marked as thick line. The 20 Myr isochrone appears to be matching the data 
reasonably well, but this pre-MS isochrone does not fully reproduce the curvatures at V — I 2. 
A similar but larger mismatch was reported by Patten &; Pavlovsky (1999; see their Figure 3). It 
seems that the Isochrones match the data better in the lower MS (cross symbols in Figure 16) 
than the fit shown in their Figure 3, which was based on a different set of isochrones. 

The level of agreement between the model and the data may not yet be satisfying to some 

readers. However, wc arc quite impressed by the matches, considering the large uncertainties still 
embedded in the convection approximation in cool stars and their color transformation. Fits to 
open clusters involve complex considerations of membership, photometry of differing quality for 
stars in differing magnitude ranges, (differential) reddening, metallicity, and of course, distance, 
all of which impinge on the age assigned. Inconsistencies are often found and resolved only in the 
context of detailed studies, and we invite the community to use these isochrones in such contexts. 

7.5. Photometric Evolution 

One of the most exciting outcomes of our isochrone tests is that the integrated colors can be 
practically used as reliable age indicators of simple stellar populations. Astronomers have long 
been aware of this fact, but the level of the accuracy of integrated colors as age indicators has 
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never been very high. Optical integrated colors were believed to be useful in determining whether 

a population is 5 Gyr or 15 Gyr, but not whether it is 10 Gyr or 13 Gyr; that it, the expected 
accuracy was of an order of 50%. This is because the color evolution of old populations is too slow 
to be useful as an age indicator. In this paper, we claim that we have achieved a substantially 
better precision. Figure 18 shows the mean integrated colors of Galactic globular clusters (Harris 
1996) compared with the models. Thin lines are the models made of only MS and RGB stars and 
thick lines are the complete models with all phases including HB stars. When complete models 
are used, the mean colors of the Milky Way clusters indicate ages of 12 - 13 Gyr, which is virtually 
identical to the age estimates obtained from the isochronc fitting. If this is not a coincidence, this 
improvement in the predictability of integrated colors as a function of age and metallicity will be 
useful to many applications in the population study. 

8. Conclusions 

If the HRD opened the window into the details of stellar evolution, it is isochrones that made 
the knowledge of stellar evolution practically meaningful to other areas of astronomy. 

In this first paper, we publish a set of isochrones in which the heavy elements are scaled to 
the solar mixture. The use of the solar mixture is quite appropriate for the study of many stellar 
populations with approximately solar relative abundances. Isochrones to the tip of the giant 

branch for ages ranging from 0.1 to 20 Gyr have been constructed for metallicities Z between 
0.00001 and 0.08. Additional isochrones with ages down to 1 Myr are also provided for the pre-MS 
and the MS phases only. 

A follow-up paper will consider mixtures that arc enhanced in a-elements as compared to the 
sun, since many pieces of evidence suggest ct-element enhancement in metal-poor stars and some 
metal-rich populations. Meanwhile, the reader is provided with a convenient correction formula 
for a-element enhancement. This formula has proved reliable in the case of metal-poor stars in 
the field and in globular clusters which are believed to be a-element enhanced. 

Important advantages of the Y'^ Isochrones over the RYI are as follows. (1) The Y'^ Isochrones 
are based on updated input physics. (2) They are available for the more up-to-date color 
transformation table of Lejeune et al. (1998), as well as for the Green et al. table used for the 
RYI. (3) The stellar models start their evolution from the pre-MS birthline instead of from ZAMS. 
This makes the Y'^ Isochrones the first comprehensive isochrones (covering a large range of age 
and metallicity) based on the stellar models that include the pre-MS stage. This fact makes the 

Isochrones useful to young open cluster studies as well. (4) Convective core overshoot has been 
taken into account. This has important impacts on isochrone fitting to open cluster HRDs and to 
the spectral dating of intermediate-age populations. 

We note that the Y'^ Isochrones match the properties of observed objects much better than 
the RYI. In particular, good matches are achieved at substantially smaller ages than before for 
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the Galactic globular clusters, by approximately 20%, once a-element enhancement corrections 
are taken into account. Approximately 15% (out of 20%) of this age reduction is caused by the 
update in the stellar models. 

Partly owing to the improved color calibration, and partly due to the inclusion of post-RGB 
evolutionary phases, the model integrated colors generated by the isochrones now also match 
within the uncertainties the observed integrated colors of Galactic globular clusters at their 
expected ages. This is one of the most significant achievements made by the Y"^ Isochrones. 

The Y"^ isochrones are available from the authors upon request or directly from our WEB 
site^, as well as electronically through this paper. We are currently constructing a new set of 
a-element enhanced Y^ Isochrones. 
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Table 1: Input Physics and Parameters 



Solar mixture 
OPAL Rosseland mean opacities 
Low temperature opacities 
Equations of state 
Energy generation rates 

Neutrino losses 
Convective core overshoot 
Alpha-element enhancement 
Helium diffusion 
Mixing length parameter 
Primordial helium abundance 
Galactic helium enrichment parameter 



Grevesse & Noels (1993) 
Rogers k, Iglesias (1995), Iglesais & Rogers (1996) 
Alexander & Ferguson (1994) 
OPAL EOS (Rogers et al. 1996) 
Bahcall & Pinsonneault (1992; 1994 priv. comm.) 
Itoh et al. (1989) 
0.2 Hp for age < 2 Gyr 
none 

Thoul et al. (1994) 
l/Hp = 1.7431 

Y = 0.23 
Ay/AZ = 2.0 
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Table 2: Chemical Compositions and Convective Core Developing Mass 



z 


Y 


[Fe/H] 


i\/fconv 


0.00001 


0.23002 


-3.29 


2.1 


0.00010 


0.23020 


-2.29 


1.6 


0.00040 


0.23080 


-1.69 


1.4 


0.00100 


0.23200 


-1.29 


1.3 


0.00400 


0.23800 


-0.68 


1.2 


0.00700 


0.24400 


-0.43 


1.2 


0.01000 


0.25000 


-0.27 


1.2 


0.02000 


0.27000 


0.05 


1.2 


0.04000 


0.31000 


0.39 


1.1 


0.06000 


0.35000 


0.60 


1.1 


0.08000 


0.39000 


0.78 


1.0 
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Table 3: Description of the Y"^ Isochroncs 



Column Description 

1 Mass in Mq 

2 logT 

3 logL/LQ 

4 log5 

5 Mv 

6 U-B 

7 B-V 

8 V -R 

9 V -I 

10 y - 

11 V-H^ 

12 V-K^ 

13 Star number density for power law IMF index x = — 1^ 

14 Star number density for power law IMF index x = 1.35'^ 

15 Star number density for power law IMF index x = 3^ 



"This color is not available in the isochrones based on the GDK table. 

''Number of stars in the box defined by this mass grid where there are 1000 stars in the IMF in the mass range of 0.5 
- 1.0 Mq. 
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Table 4: Part of a sample isochrone. The whole isochornes are available electronically. 



M/Msun 




logT 


logL /Ls 




logg 




Mv 


U-B 


B-V 


V-R 




V-I 




V-J 


V-H 


V-K 


V-L 




0.4000000 


3 


.5645 


-1.6026 


4, 


.8526 


10. 


,051 


1.076 


1, 


.506 


1.030 


1, 


.982 


3, 


,701 


3.486 


3.701 


3.701 


4 


0.4195508 


3 


.5800 


-1.5126 


4, 


.8453 


9, 


,613 


1.031 


1, 


.418 


0.934 


1, 


.790 


3, 


,429 


3.223 


3.429 


3.429 


3 


0.4391149 


3, 


.5956 


-1.4223 


4, 


.8372 


9, 


,198 


0.960 


1, 


.323 


0.860 


1, 


.656 


3, 


,136 


2.938 


3.136 


3.136 


3 


0.4586498 


3 


.6113 


-1.3322 


4, 


.8288 


8, 


,830 


0.876 


1, 


.239 


0.775 


1, 


.501 


2, 


,938 


2.754 


2.938 


2.938 


2 


0.4782554 


3 


.6269 


-1.2416 


4, 


.8188 


8, 


,492 


0.784 


1, 


.165 


0.686 


1, 


.336 


2, 


,797 


2.635 


2.797 


2.797 


2 


0.4978195 


3 


.6425 


-1.1510 


4, 


,8080 


8, 


,181 


0.719 


1, 


.102 


0.626 


1, 


.216 


2, 


,691 


2.543 


2.691 


2.691 


2 


0.5190451 


3 


.6578 


-1.0594 


4, 


.7957 


7, 


,867 


0.644 


1, 


.036 


0.570 


1, 


.103 


2, 


,570 


2.434 


2.570 


2.570 


2 


0.5409435 


3, 


.6730 


-0.9681 


4, 


.7832 


7, 


,535 


0.532 


0, 


.955 


0.507 


0, 


.982 


2, 


,389 


2.264 


2.389 


2.389 


2 


0.5633278 


3 


.6882 


-0.8773 


4, 


.7708 


7, 


,233 


0.379 


0, 


.875 


0.466 


0, 


.912 


2, 


,210 


2.095 


2.210 


2.210 


2 


0.5860503 


3 


.7033 


-0.7869 


4, 


.7580 


6, 


,960 


0.213 


0, 


.797 


0.438 


0, 


.868 


2, 


,063 


1.957 


2.063 


2.063 


2 


0.5901278 


3 


.7061 


-0.7709 


4, 


.7562 


6, 


,916 


0.184 


0, 


.782 


0.435 


0, 


.862 


2, 


,043 


1.939 


2.043 


2.043 


2 


0.5942454 


3 
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Table 5: Position of the Red Giant Branch Bump 
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Table 6: Maximum brightness (My) of the RGB as a Function of Age and Metalhcity (Z) 



t(Gyr) 0.00001 0.0001 0.0004 0.001 0.004 0.007 0.01 0.02 0.04 0.06 0.08 



2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 



-1.737 
-2.091 
-2.423 

-2.528 
-2.573 
-2.600 
-2.618 
-2.637 
-2.644 
-2.649 
-2.654 
-2.659 
-2.663 
-2.666 



-2.392 
-2.597 
-2.673 

-2.709 
-2.730 
-2.744 
-2.756 
-2.768 
-2.770 
-2.772 
-2.773 
-2.774 
-2.775 
-2.775 



-2.671 
-2.703 
-2.708 

-2.706 
-2.702 
-2.699 
-2.698 
-2.698 
-2.694 
-2.689 
-2.684 
-2.678 
-2.674 
-2.670 



-2.584 
-2.635 
-2.568 
-2.543 
-2.522 
-2.503 
-2.488 
-2.473 
-2.456 
-2.441 
-2.426 
-2.410 
-2.394 
-2.383 



-2.207 
-2.099 
-1.982 

-1.889 
-1.787 
-1.720 
-1.655 
-1.594 
-1.550 
-1.509 
-1.469 
-1.441 
-1.403 
-1.376 



-1.919 
-1.727 
-1.586 

-1.486 
-1.399 
-1.319 
-1.253 
-1.202 
-1.150 
-1.111 
-1.071 
-1.039 
-1.004 
-0.974 



-1.573 
-1.369 
-1.228 

-1.123 
-1.036 
-0.953 
-0.883 
-0.833 
-0.784 
-0.745 
-0.703 
-0.669 
-0.634 
-0.607 



-1.002 
-0.789 
-0.645 

-0.538 
-0.453 
-0.372 
-0.304 
-0.242 
-0.195 
-0.148 
-0.113 
-0.070 
-0.039 
-0.003 



-0.415 
-0.185 
-0.044 

0.083 
0.193 
0.271 
0.335 
0.393 
0.459 
0.510 
0.563 
0.597 
0.639 
0.672 



-0.052 
0.179 
0.340 

0.463 
0.564 
0.657 
0.733 
0.799 
0.839 
0.878 
0.912 
0.960 
1.007 
1.055 



0.231 
0.457 
0.622 
0.749 
0.852 
0.939 
1.015 
1.082 
1.139 
1.183 
1.217 
1.247 
1.277 
1.325 
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Table 7: Change of age estimate caused by the update in the stellar models alone 



Z Age Change (%) 
0.0001 -17.21 
0.0004 -15.87 
0.001 -17.16 
0.004 -10.57 
0.02 10.02 
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Fig. 1. — Comparison of bolometric correction (BC). The comparison is for the stars of log 5 = 4.5. 
The difference is larger in metal-poor stars. 
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Fig. 2. — Comparison in colors between the LCB table and the GDK table. The y-axis shows LCB 
color minus GDK color. Dwarfs are for logy = 4.5 and Giants are for logy = 0.5. The difference 
is substantial for giants. 




Fig. 3. — Stellar evolutionary tracks from the pre-MS birthline to the RGB, all for the solar 
composition. The pre-MS phase is shown in dotted lines. 




Fig. 4. — A set of sample isochrones for Z = 0.02, based on the LCB table . The bottom right 
panel shows the luminosity functions for all isochrones but the youngest ones of age smaller than 
0.1 Gyr, where the farthest to the right (top) is the youngest (0.1 Gyr). 
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Fig. 5. — Comparison between the Revised Yale Isochrones and the new isochrones (based on the 
LCB table) for Z = 0.02. The top left panel shows the difference in theoretical HRD which is 
caused by the difference in stellar models. The other panels show the effects of the use of the 
updated color transformation table. 
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Fig. 6. — Metal poor RGB tracks that show first dredge up bump phenomenon (see text). ^From 
bottom to top, models are for 1.2, 1.3, 1.4, 1.5, 1.6 M©. The models shown here are for no convective 
core overshoot. 
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Fig. 7. — The position of the RGB bump caused by the first dredge-up, as a function of age and 
metalhcity (based on the LCB table), (a) Each sequence is denoted by its metalhcity. Circles 
are for 2, 5, 7, 10, and 13Gyr of age from the smallest symbol to the largest, (b) Each line is 
a metalhcity sequence for given age, as denoted in the far left. Triangles are the V magnitudes 
of the RGB-Bump adopted from Fusi Pccci et al. (1990). The most metal-rich data point is for 
NGC 6553. The data for this cluster have been kindly provided by M. Zoccali prior to publication. 
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Fig. 8. — Maximum brightest in V magnitude of the RGB as a function of age and metalhcity. The 
brightest RGB means the tip of the RGB when metahicity is low {Z ^ 0.004) but a mid-point on 
the RGB when metahicity is high. The LBC table has been used for this because it reaches farther 
down to the lower effective temperature than the GDK table does. 
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Fig. 9. — The evolution of integrated luminosity for Z = 0.02 and total mass of approximately 940 
Mq in the range of 0.4 - 1.0 Mq. The smoothness of the luminosity evolution is a direct check on 
the accuracy of the mass interpolation and thus on the LF. No matter how accurately the shape of 
an isochrone may be determined, LF can easily be in error if the mass interpolation is performed 
inaccurately. Such isochrones cannot be used in the evolutionary population synthesis. 
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Fig. 10. — Stellar evolutionary tracks from the pre-MS birthline to the RGB, all for the solar 
composition. The pre-MS phase is shown in dotted lines. Filled circles are the observed MS (Gray 
1992). 
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Fig. 11. — A test to the subdwarf MS stars whose distance has been determined by HIPPARCOS 
observation. We have used the GDK-based isochrones with the following a-enhancement 
assumption: [a/Fe]=+0.4 for [Fe/H]< -1.0, and +0.2 for [Fe/H]=-0.5. 
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Fig. 12. — A test to the HRD of 47 Tuc. The HRD data have been kindly provided by A. Sarajedini. 
The estimates of [Fe/H] and reddening are from Harris (1996). [a/Fe]=0.15 adopted following the 
general pattern in the halo. 
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Fig. 13.— Same as Figure 12 but for M 68. The HRD data are from Walker (1994). The estimates 
of [Fe/H] and reddening are from Harris (1996), and [Q;/Fe]=0.3 has been adopted following the 
halo observations. 
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Fig. 14. — The use of the updated stellar models alone leads to an age reduction for the Galactic 
globular clusters by approximately 15%. A given observed MSTO brightness indicates a much 
smaller age when the Y'^ isochrones are used than when the RYI are used. The Girardi et al's 
(2000) isochrones are in good agreement with ours. 



-45- 



14 



> 




0.5 0.9 1.3 1.7 2.1 2.5 3.0 3.4 3.8 4.2 4.6 5.0 

V-l 

Fig. 15. — A test to the HRD of the metal-rich cluster NGC6553. The new isochrone (based on 
the LCB color table) matches the overall HRD well at 9Gyr. The data are from Zoccali et al. 
(2001; small dots) and from Sagar et al. (1999; crosses). The open diamond is from our synthetic 
HB models. 
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Fig. 16. — A test to the HRD of Pleiades. The data are from Schilbach et al. (1998, crosses) and 
from Eichhorn et al. (1970, diamonds). The Eichhorn et al. data are truncated at y = 10 mag, 
because they are overlapped by Schilbach et al's more recent data. Isochrones arc for Z = 0.02 and 
for ages of 1, 2, 4, 10, 20, 40, 60 (thick line), 100 (top) Myr from right to left. Adopted parameters 
are (m-M)=5.6, E(B-V)=0.04, and Ry = 3.1. 
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Fig. 17. — A test to the HRD of IC2391. This young open cluster is matched by the 20Myr 
solar-composition isochrone. The data are from Patten &: Pavlovsky (1999) and have been kindly 
provided by Patten. Adopted parameters are (m-M)=6.05, E(B-V)=0.01, and Rv = 3.1. 
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Fig. 18. — The evolution in integrated color of metal-poor populations. Compared are the observed 
integrated color of Galactic globular clusters with E{B — V) < 1.0. The data are from Harris (1996). 
Thine lines arc the models that contain only MS and RGB stars. Thick lines are the models that 
contain all evolutionary phases. The age estimates indicated by their integrated colors are in good 
agreement with isochrone-derived ages. 



